113 research outputs found
Bayesian Methods for Multi-Messenger Analysis of Supermassive Black Hole Binaries: Pulsars and Quasars and Gravitational Waves, Oh My!
Supermassive black hole binaries (SMBHBs) can lurk, often unseen, in the centers of post-merger galaxies, and pulsar timing arrays (PTAs) are rapidly approaching the sensitivities required to detect nanohertz gravitational waves (GWs) from these giant pairs. Independently, numerous electromagnetic surveys are seeking evidence of these dynamic duos’ effects on their host galaxies by searching for periodicities in time-domain observations. Combining these two methods to use multi-messenger techniques allows us to learn more about these binaries than using one messenger alone. In this thesis, we have created Bayesian methods to search for SMBHBs using electromagnetic observations of quasars and through GW emission in PTA data.
By using electromagnetic observations to identify an SMBHB candidate, we gain numerous pieces of information that also define the source\u27s GW emission, including the location of and distance to the SMBHB\u27s host galaxy and the orbital period of the SMBHB. In our study, we developed the first multi-messenger techniques used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), and applied them to a well-known supermassive black hole binary candidate, 3C 66B. We placed the lowest chirp mass limit to date on an SMBHB within 3C 66B of M
Next, we analyzed the capabilities of Bayesian methods to search for electromagnetic signatures of these binaries in simulated time-domain data sets from next-generation surveys. We developed a Bayesian model selection technique to identify periodicities induced into a quasar light curve by the orbital motion of an SMBHB from within intrinsic red noise. We discovered that future surveys, such as the Legacy Survey of Space and Time (LSST), will identify more robust SMBHB candidates than current surveys, such as the Catalina Real-time Transient Survey (CRTS).
Finally, we present the results of searches for bright continuous GWs (CWs) from individual SMBHBs in NANOGrav\u27s 12.5-year data set. A red noise process, which could be the first signs of an emerging stochastic GW background (GWB), was previously detected for the first time in this data set. In our work, we searched for CWs alongside this common noise process for the first time in real PTA data, and developed necessary data-handling techniques which will be critical for the detection of a CW, which may come soon after the potentially imminent detection of the GWB
Efficient large-scale, targeted gravitational-wave probes of supermassive black-hole binaries
Supermassive black hole binaries are promising sources of low-frequency
gravitational waves (GWs) and bright electromagnetic emission. Pulsar timing
array searches for resolved binaries are complex and computationally expensive
and so far limited to only a few sources. We present an efficient approximation
that empowers large-scale targeted multi-messenger searches by neglecting GW
signal components from the pulsar term. This Earth-term approximation provides
similar constraints on the total mass and GW frequency of the binary, yet is
times more efficient.Comment: Comments welcom
Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set
When galaxies merge, the supermassive black holes in their centers may form
binaries and, during the process of merger, emit low-frequency gravitational
radiation in the process. In this paper we consider the galaxy 3C66B, which was
used as the target of the first multi-messenger search for gravitational waves.
Due to the observed periodicities present in the photometric and astrometric
data of the source of the source, it has been theorized to contain a
supermassive black hole binary. Its apparent 1.05-year orbital period would
place the gravitational wave emission directly in the pulsar timing band. Since
the first pulsar timing array study of 3C66B, revised models of the source have
been published, and timing array sensitivities and techniques have improved
dramatically. With these advances, we further constrain the chirp mass of the
potential supermassive black hole binary in 3C66B to less than using data from the NANOGrav 11-year data set. This
upper limit provides a factor of 1.6 improvement over previous limits, and a
factor of 4.3 over the first search done. Nevertheless, the most recent orbital
model for the source is still consistent with our limit from pulsar timing
array data. In addition, we are able to quantify the improvement made by the
inclusion of source properties gleaned from electromagnetic data to `blind'
pulsar timing array searches. With these methods, it is apparent that it is not
necessary to obtain exact a priori knowledge of the period of a binary to gain
meaningful astrophysical inferences.Comment: 14 pages, 6 figures. Accepted by Ap
An unusual pulse shape change event in PSR J1713+0747 observed with the Green Bank Telescope and CHIME
The millisecond pulsar J1713+0747 underwent a sudden and significant pulse
shape change between April 16 and 17, 2021 (MJDs 59320 and 59321).
Subsequently, the pulse shape gradually recovered over the course of several
months. We report the results of continued multi-frequency radio observations
of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment
(CHIME) and the 100-meter Green Bank Telescope (GBT) in a three-year period
encompassing the shape change event, between February 2020 and February 2023.
As of February 2023, the pulse shape had returned to a state similar to that
seen before the event, but with measurable changes remaining. The amplitude of
the shape change and the accompanying TOA residuals display a strong
non-monotonic dependence on radio frequency, demonstrating that the event is
neither a glitch (the effects of which should be independent of radio
frequency, ) nor a change in dispersion measure (DM) alone (which would
produce a delay proportional to ). However, it does bear some
resemblance to the two previous "chromatic timing events" observed in
J1713+0747 (Demorest et al. 2013; Lam et al. 2016), as well as to a similar
event observed in PSR J1643-1224 in 2015 (Shannon et al. 2016).Comment: 19 pages, 8 figures. Submitted to ApJ. Data available at
https://doi.org/10.5281/zenodo.723645
Multimessenger Gravitational-wave Searches with Pulsar Timing Arrays:Application to 3C 66B Using the NANOGrav 11-year Data Set
When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational-wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C 66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C 66B to less than (1.65 ± 0.02) × 109 M o˙ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data over "blind"pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences
The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav)
has reported evidence for the presence of an isotropic nanohertz gravitational
wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by
a population of inspiraling supermassive black hole binary (SMBHB) systems,
then the background is predicted to be anisotropic, depending on the
distribution of these systems in the local Universe and the statistical
properties of the SMBHB population. In this work, we search for anisotropy in
the GWB using multiple methods and bases to describe the distribution of the
GWB power on the sky. We do not find significant evidence of anisotropy, and
place a Bayesian upper limit on the level of broadband anisotropy such
that . We also derive conservative estimates on the
anisotropy expected from a random distribution of SMBHB systems using
astrophysical simulations conditioned on the isotropic GWB inferred in the
15-yr dataset, and show that this dataset has sufficient sensitivity to probe a
large fraction of the predicted level of anisotropy. We end by highlighting the
opportunities and challenges in searching for anisotropy in pulsar timing array
data.Comment: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as
part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave
Background. For questions or comments, please email [email protected]
The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background
Recently we found compelling evidence for a gravitational wave background
with Hellings and Downs (HD) correlations in our 15-year data set. These
correlations describe gravitational waves as predicted by general relativity,
which has two transverse polarization modes. However, more general metric
theories of gravity can have additional polarization modes which produce
different interpulsar correlations. In this work we search the NANOGrav 15-year
data set for evidence of a gravitational wave background with quadrupolar
Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that
HD correlations are the best fit to the data, and no significant evidence in
favor of ST correlations. While Bayes factors show strong evidence for a
correlated signal, the data does not strongly prefer either correlation
signature, with Bayes factors when comparing HD to ST correlations,
and for HD plus ST correlations to HD correlations alone. However,
when modeled alongside HD correlations, the amplitude and spectral index
posteriors for ST correlations are uninformative, with the HD process
accounting for the vast majority of the total signal. Using the optimal
statistic, a frequentist technique that focuses on the pulsar-pair
cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6
for ST correlations when fit for separately, and median signal-to-noise-ratios
of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While
the signal-to-noise-ratios for each of the correlations are comparable, the
estimated amplitude and spectral index for HD are a significantly better fit to
the total signal, in agreement with our Bayesian analysis.Comment: 11 pages, 5 figure
The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget
Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors.
Each individual arm, composed of a millisecond pulsar, a radio telescope, and a
kiloparsecs-long path, differs in its properties but, in aggregate, can be used
to extract low-frequency gravitational wave (GW) signals. We present a noise
and sensitivity analysis to accompany the NANOGrav 15-year data release and
associated papers, along with an in-depth introduction to PTA noise models. As
a first step in our analysis, we characterize each individual pulsar data set
with three types of white noise parameters and two red noise parameters. These
parameters, along with the timing model and, particularly, a piecewise-constant
model for the time-variable dispersion measure, determine the sensitivity curve
over the low-frequency GW band we are searching. We tabulate information for
all of the pulsars in this data release and present some representative
sensitivity curves. We then combine the individual pulsar sensitivities using a
signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA
to a stochastic background of GWs, obtaining a minimum noise characteristic
strain of at 5 nHz. A power law-integrated analysis shows
rough agreement with the amplitudes recovered in NANOGrav's 15-year GW
background analysis. While our phenomenological noise model does not model all
known physical effects explicitly, it provides an accurate characterization of
the noise in the data while preserving sensitivity to multiple classes of GW
signals.Comment: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal
Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational
Wave Background. For questions or comments, please email
[email protected]
How to Detect an Astrophysical Nanohertz Gravitational-Wave Background
Analysis of pulsar timing data have provided evidence for a stochastic
gravitational wave background in the nHz frequency band. The most plausible
source of such a background is the superposition of signals from millions of
supermassive black hole binaries. The standard statistical techniques used to
search for such a background and assess its significance make several
simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often
iii) a power-law spectrum. However, a stochastic background from a finite
collection of binaries does not exactly satisfy any of these assumptions. To
understand the effect of these assumptions, we test standard analysis
techniques on a large collection of realistic simulated datasets. The dataset
length, observing schedule, and noise levels were chosen to emulate the
NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn
from models based on the Illustris cosmological hydrodynamical simulation were
added to the data. We find that the standard statistical methods perform
remarkably well on these simulated datasets, despite their fundamental
assumptions not being strictly met. They are able to achieve a confident
detection of the background. However, even for a fixed set of astrophysical
parameters, different realizations of the universe result in a large variance
in the significance and recovered parameters of the background. We also find
that the presence of loud individual binaries can bias the spectral recovery of
the background if we do not account for them.Comment: 14 pages, 8 figure
The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars
We present observations and timing analyses of 68 millisecond pulsars (MSPs)
comprising the 15-year data set of the North American Nanohertz Observatory for
Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA)
experiment that is sensitive to low-frequency gravitational waves. This is
NANOGrav's fifth public data release, including both "narrowband" and
"wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing
models. We have added 21 MSPs and extended our timing baselines by three years,
now spanning nearly 16 years for some of our sources. The data were collected
using the Arecibo Observatory, the Green Bank Telescope, and the Very Large
Array between frequencies of 327 MHz and 3 GHz, with most sources observed
approximately monthly. A number of notable methodological and procedural
changes were made compared to our previous data sets. These improve the overall
quality of the TOA data set and are part of the transition to new pulsar timing
and PTA analysis software packages. For the first time, our data products are
accompanied by a full suite of software to reproduce data reduction, analysis,
and results. Our timing models include a variety of newly detected astrometric
and binary pulsar parameters, including several significant improvements to
pulsar mass constraints. We find that the time series of 23 pulsars contain
detectable levels of red noise, 10 of which are new measurements. In this data
set, we find evidence for a stochastic gravitational-wave background.Comment: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal
Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational
Wave Background. For questions or comments, please email
[email protected]
- …